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The theoretical description of quantum phase transition, induced by the external magnetic 
field, into antiferromagnetic state in the Van Vleck - singlet - magnet with single-ion 
anisotropy of "easy-plane" type and ion spin 6" = 1 is proposed. It is shown that the 
spin polarization of the ground non-degenerated state proves to be the order parameter 
of such a transition and that the Landau thermodynamic approach can be employed for 
its (transition) description. The magnetic properties which include the field behavior of 
magnetization and magnetic susceptibility of the antiferromagnetic phase in the fields of 
different directions are studied. The peculiarities of induced magnetostriction in Van Vleck 
antiferromagnet, which as well as magnetization has a singularity in the phase transition 
point, are investigated. An attempt is made for qualitative comparison of results obtained 
with available cxpcrimantal data. 
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1 Introduction 

It is well known, that magnetization of the classical (or, what is the same, weakly anisotropic) 
antiferromagnet at low temperatures (far below Neel one Tjv), is connected with sublattice 
magnetizations turn only [1]. From this fact, it is usually supposed that their magnitude 
remains constant, and under the effect of the external magnetic field only their directions 
are changing. The character and pecuharities of magnetization process (spin-fiip, spin-fiop 
and also orientation phase transitions of the 1st kind) in such antiferromagnets depend on the 
next parameters: value and direction of the external magnetic field, anisotropy constants, 
and the magnitude of intersublattice exchange [2-5]. For example, in "easy-plane" two- 
sublattice dihalogenids NiCl2 or C0CI2 of iron group field bahavior of magnetization [6-10] 
is well satisfied by quasi-classical approach, although these magnets defer significantly. The 
value of "easy-plane" single-ion anisotropy in NiCl2 is much less than exchange - each ion 
orbital moment in crystal field is practically fully frozen. At the same time there is only 
partial freezing of orbital moment in C0CI2, and the "easy-plane" singlc-ion anisotropy is 
approximately the half of exchange (by order of its value [7]). The field behavior of induced 
magnetostriction of these crystals [11-13] are also agreed with idea of rotation of preserving 
by module of sublattice magnetizations. 

Among antiferromagnets there are, however, some family of crystalls, in which single- 
ion anisotropy exceeds inter-ion exchange [14,15] - it is so-called Van Vleck, or singlet, 
antiferromagnets. The magnetic ordering in them is absent at all temperatures, up to 
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T = 0. Such materials include, in particular, hexagonal crystals of ABX3 type, where A 
is the ion of alkah metal {A = Cs,Rb), B is the transition metal {B = Fe), X is the 
halogenyde {X = CI, Br). In these crystals magnetic moments, induced by external field 
on paramagnetic ions form antiferromagnetic chains along C3 axis, on the one hand, 

and, on the other, - triangular structures in basic plane (see reviews [16-19]). There are also 
some other compounds, that refer to Van-Vleck antiferromagnets, and among them so-called 
DTN, which chemical formula is NiCl2-4SC(NH2)2 [20-23]. It also has (antiferromagnetic) 
chains Ni-Cl-Cl-Ni along "heavy" magnetic axis, although in field absence the mean spin 
on each site is equal to zero, because of exceeding by single-ion anisotropy of parameters of 
both intra- and intersublattice exchange. It should be emphasized, that DTN can be refered 
to the group of two-sublatticc Van Vlcck antiferromagnets, that have different from NiCl2 
crystal structure, another character of exchange interactions, which by the value are much 
less than single-ion anisotropy [22,23]. 

The magnetization process in Van Vleck antiferromagnets differs principally from that 
which takes place in classical Neel antiferromagnets [25-28]. At first, the magnetic ordering 
without external magnetic field is absent and therefore there are no any magnetic sublattices. 
Secondly, magnetic, or particularly antiferromagnetic, ordering in Van Vleck antiferromagnets 
can appear spontaneously by quantum (in definition of Ref. [24]) phase transition, induced 
by magnetic field [16-23]. 

Weak dependence of magnetic susceptibility on external magnetic field, thus, represents 
some sort of peculiarity of anriferromagnetic phase. As a result, the observed magnetization 
actually follows the hnear field behavior [22,23,29,30]. It appears, in other words, that such 
a behavior of magnetization of a system (but not the proper sublattice magnetization) in 
antiferromagnetic phase is similar to the magnetization induced by external field in Neel 
antiferromagnets. It could be understood, if the transition to the antiferromagnetic phase 
would be the phase transition of the 1st kind. In such a case sublattices could magnetize, in 
the transition point, due to jump (in the presence of corresponding susceptibility singularity), 
and at further field growth there can occur sublattice magnetic vectors turn only. The 
experiment shows, however, that transformation of non-magnetic (singlet) state in the 
antiferromagnetic one takes place continuously, what means, that this magnetic transformation 
is the phase transition of Ilnd kind [22,23,29,30]. The latter demonstrates, that sublattice 
magnetizations appear and change also continuously from their initial zero value to maximal 
one. So, the classical approach with the constant module of sublattice average spin for Van 
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Vleck antiferromagnets is not applicable fundamentally. 

The results of DTN induced magnetostriction measurements are presented in the papers 
[22,23]. It could be seen from them, particularly, that induced magnetostriction in DTN 
appears and exists namely in the antiferromagnetic phase only. There was also found, that 
relative crystal deformation, that arised along "heavy" magnetic axis, changes its sign 
during field increasing. Such a behavior of DTN striction was attributed in Ref. [23] to the 
prevailing role in this compound of intersublattice magnetoelastic interaction of exchange 
nature [20,21]. On the other hand, the magnetostriction sign change is also observed in some 
classical Neel antiferromagnets, for example in C0CI2 [11-13], where it is conditioned by the 
anisotropic intrasublattice magnetoelastic interaction. 

Thus, the description of induced magnetostriction in the Van Vlcck antiferromagnets, 
where anisotropy is not small, and so, anisotropic magnetoelastic interaction also can be 
comparable (or even larger) with isotropic exchange one, becomes relevant. The consideration 
of this problem requires to account the fact, that sublattice magnetization modules in 
antiferromagnetic phase of the initially singlet magnet essentially depend upon the field. 

From the above said it can be evident, that there are some unresolved questions of the 
description of induced magnetic phase transition into the antiferromagnetic phase and its 
magnetic characteristics (field dependencies of sublattice and system as a whole magnetizations, 
magnetic susceptibility, magnetostriction) in Van Vleck systems. 

Below we shall proceed from the suggestion, that in Van Vleck magnets the intrinsic 
spontaneous magnetic (or antiferromagnetic) moment is equal to zero, so for them there is 
no magnetic ordering temperature without the external magnetic field. It would seem, that 
the absence of magnetic (dipole) moment , or, in other words, magnetization (spin) on the 
site, shows not only the absence of any magnetic ordering, but, clearly, the absence of the 
magnetic contributions in physical properties of corresponding systems. However, in fact, 
this is not right, because the absence of ordinary - exchange- induced - spin ordering do not 
include the presence the ordering of other type in, particularly, the quadrupole one. The 
latter, in one or other way, is peculiar to all Van Vleck magnets, which are the special case 
of magnetic crystals with more specific - nematic - type of spin ordering [31,32]. 

The one or another spin ordering proves itself not only by the appearance in crystal of 
new (spin-)electron excitation branches, which, for example, in Van Vleck nematics turn out 
gapped. It can also be revealed in such an observed and calculated characteristic of these 
magnets as their magnetostriction, which peculiarities for such systems is not completely 
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studied yet. A the same time, because of recent calculations of DTN magnetoelastic features 
[22,23], there is definitely such a need. 

It can also be noted, that in some papers (for example, Refs. [33-36]), the description 
of phase transition between singlet and induced atiferromagnetic states carries out by using 
the representation of Bose-Einstein condensation of magnons. Indeed, the appearance of 
magnetization in finite magnetic fields can be formally described in the terms of some 
magnetic excitations condensation. But in reality in observed systems there does not occur 
any true condensation of quasiparticlos, because, as it will further be seen, one should say 
about rearrangement of the ground state only, and thus - about virtual, but not real magnons 
[37]. 

Below there is considered the model of strongh^ anisotropic, two-sublatticc antifcrromagnct 
with ion bare spin S = 1. In the framework of quantum approach it will be made an 
attempt to describe the crystal magnetization, magnetic susceptibility and magnetostriction 
at magnetically induced phase transition from the initial singlet state to the spin-ordered 
one. For calculation of physical characteristics of a system, there will be used the total energy 
E, that is the sum of relevant contributions: 

E = Eexch + Ean + Eh + E^l + Em-el, (1) 

where Eexch is the exchange energy; Ean is the magnetic anisotropy energy; E^ is Zeeman 
energy; Eei is the elastic energy and E,n^,,t is the magnetoelastic energy, or the energy of 
spin-lattice coupling. It is also supposed, that magnetoelastic interactions arc much weaker 
then exchange ones, and do not have a noticeable feed-back influence on the magnetic 
ordering. Assumptions made allow to provide a calculation in the simplest, but rather general 
form, confining, as usual, by summands, that are linear by elastic deformation tensor in 
the magnetoelastic energy and are quadratic by this tensor in elastic energy. So, at this 
hmitations, there can be solved a problem of the magnetic ordering, at first, and only then 
use obtained results for field dependence of induced striction. 
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2 The ground state of singlet antiferromagnet with S = 1 
and "easy-plane" magnetic anisotropy in the longitudinal 
magnetic field 

In accordance with above mentioned the consideration will be for simplicity restricted 
by bihneal anisotropic (intra- and intersublattice) exchange interactions, single-ion "easy- 
plane" anisotropy and Zeeman term. The simplest model Hamiltonian of a system, that 
defines three contributions, Eexch, Ean and E^, in Eq. (1), can be written as: 

where a, /3 = 1,2 are the magnetic sublattice indices, which numbers in the considered 
system was chosen as 2; vectors n and m gives spins position in magnetic sublattices, which 
are described by spin operators Sno,; the constant D > 0, that reflects an "easy-plane" 
magnetic structure; the magnetic field /;, is defined in units of energy, so h\\ = HBfjHz] Hz 
is the OZ projection of magnetic field, at that the crystallographic symmetry axis OZ is 
perpendicular to the "easy" plane. Just at the h || OZ the magnetic field induce the phase 
transition to the antiferromagnetic state. The case of transverse fleld h _L OZ will also be 
considered below, but it should be emphasized, that such a field do not induce any phase 
transitions. Parameters Jncm^ characterize the value of exchange interaction isotropic part 
and Ja„mi3 is exchange anisotropy, which can be basically both "easy-axis" and "easy- 
plane" types. We, however, will assume, that inter-ion anisotropy, as also single-ion one, 
relates to the same - "easy-plane" - type. 

The convenience of these restrictions is conditioned by the fact, that in such a physical 
situation both sublattices become symmetric relatively to the external field, what allows to 
reduce twice the number of equation derived. 

The analysis of possible quantum eigenstates of Hamiltonian (2) at h || OZ will be 
provided, using the approximation of self-consistent field, that corresponds to spin fiuctuation 
neglecting and to the change of average by multiplying spin operators of different sites on 
multiplying of averages. In this case the energy Egr of the ground state, normalizing on one 
cell (for nearest both inter- and intrasublattice different spins) is equal to: 



a/3 a/3 a a 

where Sa is the quantum-mechanical average of spin vector of ath sublattice in the ground 
ion state; Zaf} is the number of nearest neighbors from the same {zaa) and another {Zaf} = 
Z12) sublattices. Also here are introduced such averages for components of spin quadrupole 
moment Qf^ [40-42]. It should be also noted that for antiferromagnet the intersublattice 
exchange parameter is J12Z12 = / > 0; at the same time, the parameter JnZn = J22Z22 = 
J of intrasublatticc exchange can be of any sign, that we will also choose for simplicity 
as furthering to the ordering, J < 0. The exchange anisotropy, in this case, satisfies the 
conditions of its "easy-plane" type: Ji2Z\2 = AJ < and Jn^ii = J^i^i^ = A J > 0. 

Let us impose for spins of each sublattice their proper (rotating) coordinate systems 
Ca^aCa) that ath sublattice average spin is always oriented along axis, what means that 
this axis is the quantization one for this spin sublattice, and axis is lain in ZC^a. plane. 
Thus, in such coordinate systems the correct wave function of the ath sublattice ground 
spin state, as well known, has next form [38,39]: 

■0^°^ = cos?!>a I 1 > +sin^„ I -1 >, (4) 

where | ±1 > and | > are cigcnfunctions of operator 5^ ^ in bra-kct representation. Next, 
it can be calculated, using (4) the quantum-mechanical spin and quadrupole averages: 



s = cos20, 0" = 1, = - (H-sin2(?;.). (5) 

In the expressions (5) the sublattice indices are missed, because as mentioned at the 
chosen field direction the evident dependence of observables on index a is absent. 
The usage of functions (4) allows one to describe the energy (3) at h || OZ as: 



Egr = I cos 2(/) cos 26* - I J| cos 2^ - Jz cos 2^ cos Q 
,2 a 



+2D 



cos^ e + 



sm 



■ (1 + sin 20) 



2h\\ cos^cos2(^. 



(6) 



where Jz ^ AJ ~ AI. 

For the determination of magnetization, magnetic susceptibility and subsequently striction, 
there should be found field behavior of mean spin (and its direction) for each sublattice in 



the field. Also it should be made the same calculation for spin quadrupole moment. As it was 
reported in Refs. [40,41], the solution of the problem of spin configuration in the magnetic 
field consists of the minimization of the expression (6) by all available unknown quantities: 
the geometrical angle 6 and (see Eq. (4)) the angle (j) of quantum states mixture. Such a 
method of observables finding, being completely an equivalent to the solution of quantum 
self-consistent problem, is more convenient and consistent, because allows the generalization 
on the case of finite temperatures [27,28]. 

The equations for both required angles are: 



del) 



= -2 {lcos29- \J\ + Jzcos^e) sin40 + 2£)sin2 6'cos20 + 4/i||Cos6lsin2^!) = 0, (7) 
dE 

-—^ = -{21+ Jz) sin 29 cos^ 2(1) - D sin 29 (1 - sin 20) + 2/i|| sin 9 cos 20 = 0. (8) 

It is known from Ref. [37], that the set of Eqs. (7)-(8) in the absence of the external 
magnetic field has two solutions: the non-magnetic one, s = 0, that exists at £> > 2(/ + | Jj) 
and the "magnetic" one at D < 2(7+ |J|), with which the reduced value of single-site mean 
spin 



s=Jl 2<1 (9) 

V + 

is associated. 

The initial ground state of the system should be the singlet one, s = 0, that the quantum 
phase transition (at the magnetic field h || OZ) from this state to magnetically ordered one 
occurs. So, below we will suppose, that the next evident inequality 2{I + | J|)/D < 1 takes 
place. At this model parameters ratio, the ground state of the system is really nonmagnetic 
and the ordering in the absence of magnetic field is not realized at any temperatures [41]. 
Otherwise, this ratio actually gives the condition on singletness of magnet ground state, 
which is Van Vleck one. The solution s = also satisfies Eqs. (7)-(8) for some interval of 
magnetic fields. 

As field grows the finite value, s 7^ 0, of the mean spin on the site appears. It can 
be derived from Eq. (8) the expression for the average spin orientation relatively to the 
crystallographic axis: 

cos 20 

^ " D (1 - sin 20) + (2/ + Jz) cos2 20 ' 
It is seen from Eqs. (7) and (10), that in large fields, when ^ hfup (where /i/Hp = 
D + 21 + J z) the state, in which spins of both sublattices are directed along "heavy" {9 = 0) 
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axis, is established. Then the spin projecton on the external field direction will be maximum 
and equal to s = S = 1. For < hfUp spins of sublattices are orientated at finite angle 
< ^ < '7r/2 to the "heavy" axis. 

3 Thermodynamic analysis of spin states 

Using the formulas (5) and substituting Eq. (10) into the Eq. (6), it could be obtained the 
ground state energy in the form of functional 

Egr = -(!+ \J\) s^ + D(l-y/l-sA , \ ^ , (11) 

which depends on the spin polarization s only. The expansion of this energy over the small 
s gives: 

^■^ li u \ 2 D 2hl{2I+Jz)\ 4 

E,r = JJ {hs -H)s' + jil+ -^^^ J (12) 



where hg = D^yl — 2{I + \ J\)/D is the critical field of magnetization appearance. 

In the expansion (12), which refers to the field region — > hg, one can restrict by 
terms, that are not higher then of 4th power by s. Actually this expansion for the ground 
state energy is similar to the free energy expression in Landau theory of phase transitions. 
However, in Eq. (12) the ground state spin polarization corresponds to the order parameter, 
and the leading parameter, that results in the phase transition, is not the temperature, but 
the magnetic field. It can be seen from Eq. (12), that at /ly < hs the coefScients at and 

are positive, and so the ground state of spin system will be Van Vleck non-magnetic 
single-ion state. At the pint = hs the sign of coefficient at changes, and in the fields 

> hs the spin polarization (of still non-degenerate ground state) spontaneously appears. 
The value of polarization can be readily found by minimization of Egr (12): 



ds ~ ' 



D ^ " 4 V 



= 0. (13) 



From Eq. (13) it follows, that near the critical point > hg this polarization (or simply 
the spin of the ground state) fundamentally depends on the field: 



/ ^ / (hu - hs) 

''W'\/ D^ + 2/.;(27 + i)/D - 

In the same vicinity, /;,|| > hs, of the critical point the angle between vector s and axis 

OZ is determined by the expression: 



S 



COS0 



2D\ D-^ + 2Ki {2I+Jz)/D' 



(15) 



Thus, it is found, at h\\ ~ hg the spin polarization spontaneously arises as field grows 
in the very "easy" plane, because at — hs ^ the angle 6 7r/2. In other words, it 
turns out, that in the moment of its appearance, the vector s(/i|| > hg) is perpendicular to 
the longitudinal field: s ± H || OZ. Further magnetic field growth leads not only to the 
decreasing of 0, at it follows from Eqs. (14)-(15), but also to the simultaneous increasing of 
spin polarization, that is as bigger its value, as more it flatten against the "heavy" axis. 

In the whole, the induced tilt of the magnetic subalttices, and thereafter the magnetization 
of Van Vleck antiferromagnet include two processes: the classical rotation of spins (sublattice 
magnetizations) and purely quantum (because of angle (j) change) growth of single-site 
polarization s(/i||). Both processes take place also at T = 0. The antiferromagnet magnetization 
(normahzing on one magnetic atom) is described by the evident product: 



As a result, one arrives to the unexpected result: the observed magnetization near 
critical fleld of quantum transition from singlet to spin-polarized state depends linearally 
- as in classical antiferromagnets - upon the external magnetic field, that induce the very 
transition. From this comes another rather remarkable conclusion: at such a phase transition 
the magnetic susceptibility of a system should have a jump. 

Thus, in the framework of approach, that is similar to the Landau thermodynamic one, 
it was demonstrated, that the spin polarization is the only order parameter for induced by 
magnetic field h || OZ quantum phase transition from Van Vleck phase to the antiferromagnetic 
one. But despite the fact, that calculation was made for the case T = 0, the required 
polarization proves to be essentiahy dependent on the external field. It should be reminded, 
that in classical antiferromagnets ion spin polarization is fixed at T = and do not evaluate 
in the field, while in Van Vleck system it appears as a consequence (in terminology of Ref. 
[24]) of quantum phase transition [27,28]. 

Next an attention should be drowned to such an analogy, that single-ion anisotropy, 
reducing average spin, plays a role of "disordering" factor, and in this sense can be 
compared with entropy. It (single-ion anisotropy) leads to the mixture (or linear combination 
of quantum states), that results in the absence of spin polarization of ions in their ground 
state. The exchange and magnetic fields, on the contrary, resist to this, "magnetizing" the 





(16) 



D3 + 2/i2 {21 + Jz) 



system and causing the spontaneous (or forced) spin polarization, which at the moment of 
its appearance is directed perpendicularly to the magnetic field. 

The studied quantum phase transition between Van Vleck (also ordered, in essential) 
and antiferromagnetic states is, as it was seen, the consequence of different interactions 
(exchange, Zeeman and spin-orbital, that lies at the heart of single-ion anisotropy) competition. 
Therefore such a quantum transformation is natural to identify as the magnetic phase 
transition of "displacement" but not of "order-disorder" type. As distinct from the last 
one, the transition of displacement type is not the transition in the system of spins, which 
fluctuate between degenerated (or almost degenerated) quantum states "up" and "down because 
the ground state of quantum paramagnets is always non-degenerated and its polarization is 
the direct consequence of this state rearrangement in the external field. 

There should be noted, at last, that the apphcability of phenomenological theory, that 
is based on the expansion (11), is confined by the fields /i|| > hs in the vicinity of critical 
point hs- In the field region » hg the magnetization process should be analyzed with 
more exact expressions both for ground state energy and for the equations, which define the 
spin configurations. However the latter can be easily found numerically. 

4 The magnetization curves and magnetic susceptibility 
of antiferromagnetic phase 

Substituting in the Eq. (11) the expression (5) the equation, that describes the spin polarization 
as the function of longitudinal field, can be obtained: 

( I ^ D{\ + ^Y^)h\ \ 

s £»-2(7+|J|)Vl-s2 —= — 2 =0- (17) 

V (iP(l + x/r^)+(2/ + Jz))V 

It should be noted, that this equation refers both the fields /ly < hg of the existence of 
Van Vleck phase, where (see Eq. (12)) = ■^D'^ - 2£)(7-F |J|), the nonmagnetic, s = 0, 
state is stable, and to the region hs < < hfUp of the antiferromagnetic phase (up till the 
point hfUp of its fiipping). It is obviously, that at the point (see Eq. (10)) /ly = h/up, which 
is corresponded to the value ^ = 0, the polarization arrives to its maximum value s = 1 on 
the site. 

Using Eq. (17), the behavior s(ft,||) in the fields hs < < hjnp can be found, and from 
Eq. (10) - also angle 9. After that, it is not difficult to define the field dependence of the 
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Phc. 1: Magnetization mjj (solid curves), spin s, quadrupole moment Q and value of cos6' 
versus field. The curves for s, Q and cosd are calculated at \ J\/D — 0.455 and at condition 
I = Jz = 0. Magnetization lines 1-5 are for next parameters: line 1 is for \J\/D — 0.455, 
I ^ Jz = 0; line 2 is for \J\/D = 0.35, I ^ Jz = 0; line 3 is for \J\/D = 0.35, / = 0, 
Jz/D = 0.5; line 4 is for \J\/D = 0, I/D = 0.35, Jz = 0; line 5 is for \J\/D = 0.15, 
I/D = 0.2, Jz/D = 0.5. 

quadrupole ^ in Eq. (3) . In the framework of such an approach the field dependencies of 
magnetization were calculated (see Fig.l). 

The curve 1 on Fig.l was plotted for such case, when intrasublattice exchange is prevailing, 
while intersubalttice and anisotropic ones are omitted. At chosen parameters the magnetic 
sublattices are fixed artificially, because this extreme case corresponds actually to the two 
independent antiferromagnets. The field dependencies for s, Q'^^ = Q and cos^ are also 
shown on Fig.l for these parameters. . It can be seen, that at the point = hg the average 
site spin really spontaneously appears, and exists in the region > hs. With further field 
growth, the value of s is increasing, leading by velocity of the change of angle 9. This 
velocity, however, becomes more fast, when the field approaches to the flipping fleld and, 
correspondingly, s — > 1. In such a case the intrasublattice exchange (due to its isotropy) has 
no effect on spin saturation, and so the value hfup of this critical field is completely defined 
by the single-ion anisotropy. 

From curve 1, that refers to the case D — 2\J\ << D, it follows, that the fields hfUp 
and hg differ quite weakly {hfup/hg ~ 3), although in the experiment [22,23] their ratio 
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reaches 6, and from the data of Ref. [30] this ratio is about 8. Besides, the field dependence 
of magnetization for the considered case 7 = reveals, as can be seen, large nonlinearity, 
while the experimental data for all above mentioned compounds, are evidence of near linear 
dependence of magnetization on magnetizing force. 

It should be noted, that the case, in which inequality {D — 2\J\)/D « 1 is satisfied, 
is physically available, but it can not be justified from the experimental point of view. To 
demonstrate this on Fig. 1 the curve 2 is plotted, for which the difference between parameters 
in intrasublattice exchange and single- ion anisotropy is chosen not less, but bigger than for 
curve 1. This choice really leads to the increase of the field value hs and do the decrease 
of the ratio hfup/hs, what indicates, that in attempts of interpretation of the experimantal 
magnetization, the intcrsublatticc exchange can not be neglected. 

It is interesting, that when I/D 0, Eq. (17) has an exact solution, using which the 
ground state energy can be written in the form of function of magnetic field: 

Then the magnetization (normalizing on one magnetic ion again) takes next form: 

The dependence (19) is described by the lines 1 and 2 on Fig. 1. From Eqs. (18) and (19) 
it can be also seen quite a big field nonlinearity of magnetization in the antiferromagnetic 
phase. The expression (19), for fields hg can be also presented in the form of Eq. (16), 

when 2\J\ D. 

Now let consider the opposite limiting case, when intersublattice exchange is the biggest 
one in the system. It is seen, that even if one preserves the exchange field (which formahy 
gives the same value of hs), which affects on the spin from other sublattice, the change of 
magnetization (curve 4 on Fig. 1) occurs. The antiferromagnetic (intersublattice) exchange, 
unhke the intrasublattice one, leads to the growth of the field hfUp, because in this case the 
external field should overcome the effect of the same anisotropy, from one side, and, from 
the other side, - of exchange field, that prevents parallel orientation of sublattice spins. 

The curve 3 already shows the nonhnearity decrease in m{h\\), which as if it is rectified by 
intersublattice exchange (or by its anisotropic part). At the same time, the antiferromagnetic 
exchange together with the external magnetic field (in the region h\\ > hg), resisting the 
anisotropy, leads to the estabhshment of spontaneous polarization. In the large fields, when 
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polarization tends to its maximum value, the behavior of exchange in Van Vleck antiferromagnet 
does not differ from that one in classical antiferromagnets: it simply resists to the parallel 
configuration of both sublattice spins. 

Curves 3 and 5 demonstrate the influence of easy-plane exchange anisotropy. Actually 
this anisotropy does not change the position of critical fleld hg, but it also does not "want" 
the establishment of collinear state, when Si — > S2 || OZ. At the same time the account of 
exchange anisotropy of easy-plane type allows to obtain such a behavior of magnetization, 
that is near to linear one and observed in Refs. [22,23,29,30]. For clarifying how good the 
hnear dependence corresponds to m{h\\), there is shown on Fig. 2 the magnetic susceptibility 
X\\ = dm,\\/dh for the same parameters as on the Fig. 1. 

Because the magnetization is nothing other, then m[| = scoad, where s = s(/?,||), the 
longitudinal magnetic susceptibility of Van Vleck magnets is naturally to represent in the 
form of two above mentioned terms - the classical Xci and the quantum Xguan ones, so that 

X|| ~ Xcl + Xquan't 

■ „ds , , 

Xc/=SSm6^— ; Xquan=COS0—. (20) 

As can be seen from Fig. 1, near hs the biggest growth reveals the spin polarization 
s, so in the flelds /i|| —^ hg the "quantum" contribution will dominate in the magnetic 
susceptibility. And when the value of spin polarization saturates {s{h\\ — > /i/Hp) — * 1), the 
susceptibility will be mainly controlled by classical term (see Eq. (20)). 

Figure 2 shows, that when intrasublattice exchange is really largest one, then magnetic 
susceptibility grows, changing in fleld region hg < /i|| < hfUp in four times. If intersublattice 
exchange and/or exchange anisotropy "switches" on, then fleld dependence of differential 
magnetic susceptibility X|| = x(^||) becomes essentially weaker. Nevertheless, it can be seen 
from plots, shown on Fig. 1, that nonlinearity of function m{h\\) in the antiferromagnetic 
phase at the chosen parameters remains quite noticeable. 

The case, when within the boundaries of this phase the value oix{h\\) changes (30-50%), 
is shown on Fig. 3, which meets the model parameters \J\/D = 0.05, I /D = 0.3, Jz/D ~ 1 
or Jz/D = 1.5. In other words, the exchange anisotropy is comparable or even exceeds the 
single- ion one. At such ratios between the parameters the field hfHp is almost in five times 
exceeds the fleld hs (one should note, that experimentally observed ratio is h/up/hs ~ 6 
[22,23]). 

On the same Fig. 3 the dependencies of s(/i||), cos^(/i||) and (5(/j||) are shown for 
parameters \J\/D = 0.05, 1/D = 0.3 and Jz/D = 1. The behavior of (5(/i||) almost coincides 
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Phc. 2: Magnetic susceptibility x\\ versus field. The numbers of curves are corresponded with 
parameters, that were used for the plotting of lines with the same numbers, that are on Fig. 
1 

(see Fig. 1) with the field dependence m(/i||). Moreover, it follows from Fig. 3, that exchange 
anisotropy, even a comparable with single-ion one, do not fufiy linearize the function m{h\\). 
As was mentioned, this fact can be explained by the existence of two different regions in the 
magnetization of Van Vleck magnet. 

On the first of these regions, near hg, the quantum process, as was pointed out, is 
determinant and the magnetization is defined basically by the appearance and growth of 
s(/i||). On the second one, in the vicinity of /ly < hfup, the more important becomes the 
classical rotation of sublattice spins to the field direction, at essentially less (but not absent) 
role the very spin value change. It is obvious, that in this region the susceptibility depends 
much weaker on the value of magnetic field. So, it can be supposed, that the fiipping of tilted 
antiferromagnetic sublattices, or the transition of Van Vleck system between induced two- 
and one-sublattice magnetically ordered states occurs as an orientation phase transition in 
ordinary antiferromagnet, when the variation of sublattice magnetization directions is in 
fact the only process. 

However, even for this field transition the quasiclassical approach do not give the correct 
solution for m{h\\) in spin nematic. Indeed, one could suppose, that near the fiipping field, 
when s(/i|| — > hfUp) w 1, the quasiclassical magnetic energy in the ground state takes the 
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Phc. 3: Longitudinal magnetization m|| and magnetic susceptibility x\\ versus field at 
\J\/D = 0.05, I/D = 0.3, Jz/D ^ I and 1.5. The functions s{h\\), cose{h\\) and Q{h\\) 
are sliown only for | J|/D = 0.05, I/D = 0.3, Jz/D = 1. 

form: 

Egr = I COS 20 + {Jz + D) cos^ e - 2/i|| cos 6* = 0. (21) 
Then from Eq. (21) immediately follows the equation 

dE 

— ^ = 2\- (27 +Jz + D) cose + /ill] sin 61 = 0, (22) 
a/i|| 

which shows, that in the vicinity of h\\ hjup the magnetization of tilted {0 ^ 0) phase is 
proportional to the field: m|| = X||/i||, where 

= D + 27+J. " 7^ = 

It can be scon, that in the field /i/;ip, the magnetization (on one spin) is = 1. However, 
the susceptibility (23) differs from the exact ratio (20) and gives physically incorrect behavior 
of magnetization. It is connected with the fact, that in the region h\\ h/up it appears, 
that my depends linearly versus magnetic field, and asymptotically tends to zero at /ly 0. 
As for plots, shown on Figs. 1 and 3, it is easy to see, that the function m(ft,||), although 
it behaves linearly by field, but nevertheless it depends on magnetic field in not a direct 
proportion. 
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An important conclusion follows from this: the quasiclassical approach (21), based on the 
substitution of quadrupole moment Q^^ by the average spin Zth projection square, appears 
to be unapphcable even in such a field region, where spin polarization almost reaches its 
saturation value s — > 1. 



5 Field behavior of magnetization in transversal field 

As it was reminded, the phase transition to the antiferromagnetic state does not occur at 
h _L OZ, although the magnetic field magnetizes the system. 

Lets suppose, that due to the antiferromagnetic exchange in the easy plane, two sublattices 
are formed. Then their spins lay in this plane and are identically tilted to the field. In this 
case the energy of the ground state is: 



Egr = I cos 2ip cos^ 2(j)-\J\ cos^ 2<j) + D{l+sm 2(f)) - 2h± cos if cos 2(f), (24) 

where (p is the angle between vector Si (or vector S2) and field h, and the angle between Si 
and S2 is twice larger, 2(p. 

The spin configuration will be defined, as always, by minimizing the energy (24). As a 
result, it is the set of equations (cp. Eqs. (7) and (8)): 



— = - 27 cos^ 2(f) sin 2ip + 2h± sin cos 2(5!). (25) 
dip 

dE 

= -2 (/cos 20- |J|)sin4(/) + 2£)cos20 + 4/i_LCOS(/?sin2(/) = (26) 

The equation (25) has two solutions. For the first of them <^ = and it corresponds to 
one-sublattice magnetization, when the polarization of magnetic ions is directed along the 
field. The second solution cosyi = /i_l/(2/cos2(^) provides the existence of two sublattices. 
The last, taking into account Eq. (5), can be rewritten in the usual form: 



cos(^ = /i_l/27s. (27) 

The denominator of Eq. (27) is the intersublattice exchange field, and this expression is 
similar to the expression for the magnetic sublattice tilt angle in classical antiferromagnets 
[1,2]. Nevertheless, it should be noted, that spin in Eq. (27) is not equal to its maximum 
value. 
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Phc. 4: Field behavior of m{h±) (lines 1-4), Q^^{h±) (line 5) and x± for h _L OZ. The line 
1 is calculated for parameters \ J\/D = 0.455, / = 0, the line 2 is for \ J\/D = 0.35, / = 0, 
the line 3 is for \ J\/D = 0,1 = 0.35, lines 4-6 are for \ J\/D = 0.05, / = 0.3 

Substituting (27) in (26), one arrives to the equation: 

2[2(/-|-|J|)sin2(/)-hD]cos2(/) = (28) 

It follows from Eq. (28), that the spin polarization for antiferromagnetic state at h _L OZ 
should be equal to (9). However, the model parameters, accepted above, are such, that the 
denominator under the root in Eq. (9) is larger than 1, and non-polarized singlet is the 
ground state of ions. Thus, the solution (27) is possible only for initially polarized ground 
state, or when the antiferromagnetic (not singlet) phase is realized in the system even at 
h± = 0. But if without field the polarization is s = 0, then from the set of Eqs. (25)-(26) 
fundamentally another result follows: the critical field of polarization appearance in the 
transversal geometry is the h± = 0. The distinction from "longitudinal" case, for which the 
critical field is finite, is easy to explain. At any fields h± ^ the ground state with Sz = 
(in crystal coordinate system) is immediately admixed with the ionic exited states, which 
have Sz = ±1- In the case of longitudinal field, there is a threshold for such an admixture. 
Then taking into account, that the transversal field does not induce the antiferromagnetic 
phase, one can obtain: 
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Egr = {I- I J|) + £> (l - Vl - s^) - 2/i±s. (29) 

It should be also noted, that at h ± OZ the spin polarization is always equal to magnetization, 
which is directed along h, i.e. m± = s. Minimizing energy (29), one readily arrives to the 
equation 

^^'"^ =2{I-\J\)s + D^J==-2h±=0, (30) 



ds ' ' " 



s 



that allows to define the dependence of spin polarization upon the transversal field. 

On Fig. 4 there is shown the field behavior for m±, that are obtained from Eq. (30). It 
is seen, that if intrasublattice exchange dominates in the system, then the magnetization 
increases rapidly, and if intersublattice exchange is "added then the magnetization slows 
down. 

Despite the fact, that at h _L OZ the average spins are oriented also perpendicularly to 
Z, the spin quadrupole moment Q^^ versus field reveals the behavior (sec Fig. 4) similar to 
the spin polarization. The magnetization rate decreases as the field grows and the magnetic 
susceptibility has a maximum at h± 0. The normalized magnetic susceptibility x^''' = 
x{h± — 0) is also plotted on Fig. 4. 

Using Eq. (30), the expression for magnetic susceptibility at h _L OZ can be obtained; 
it has the form: 

X± = X {h±) = (31) 

2(/- |J|) + D(l-s2)-3/2 

It is seen, that in large fields, when s — > 1, transversal susceptibility x± ~* 0- In the opposite 

limit h± — > 0, the magnetic susceptibility is equal to: 

= D + 2{I-\J\) ^^^^ 
It also follows from Eq. (32) , that when the intrasublattice exchange increases the value 
of x^±^ grows and, on the contrary, at the increasing of intersublattice exchange it decreases. 
It is usual situation for physics of phase transitions, because the growth of intrasublattice 
exchange at / = can result in ferromagnetic state with characteristic to such kind of 
transition the susceptibility singularity (it goes to the infinity at the transition point). At 
the same time, the transition to the antiferromagnetic state is not accompanied by the 
abnormal growth of magnetic susceptibility. Indeed, the point of phase transition to the 
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antiferromagnetic phase corresponds the equahty D = 2{I + \J\). At its substitution in Eq. 
(32) the value x±^ = 1/2/ is directly obtained. The same will be the value of magnetic 
susceptibility in the antiferromagnetic phase, the magnetization of which is determined by 
the expression (27). 

6 Induced magnetostriction in the longitudinal magnetic 
field 

Considering the striction properties of Van Vleck antiferromagnets, it will be for certainty 
supposed, that the crystal has a hexagonal structure. There will be confined, at that, in 
magneto-elastic energy (see Eq. (1)) by spin-deformation interaction, that is proportional 
to second power of average spin [45]. Besides the single- ion terms will be accounted in the 
energy Em-el', they contain the average components of spin quadrupole moment [46,47]. 
Then for such approaches the elastic and magneto-elastic contributions to the total energy 
(1) can be presented as following: 

Eel = ^Cii {ul^ + uly) + ^CS3UI^ + CuUxxUyy + C13 {u^^ + Uyy) U^^ 

+2cii {ul^ + ul^) + 2ce6uly (33) 

Em-el = ^ [KffUzz + 7a/3 i^xx + Uyy)] SaSff + ^ [^n"'^ {Qa'^Uxx + Qa^Uyy) 

a(i a 

+AB^^r^Q^^u.y] + Y.{B[f^ {s^sfu^x + sls^.uyy) + Bf^s'^sluzz (34) 

) {s^sjuyy + slslu^^)+4.B^^^\slsfuyz + s^s$u^z)+4.B^^^\^s}u^y], 

where XajS, 7a/3 are the parameters of magneto-elastic exchange interactions, in which indices 
a, (3, as above, are the numbers of the spin sublattices, -B^p*^ and B^^^^ are the parameters of 
anisotropic magneto-elastic interactions [45], where the upper index shows on the single-ion 
or ionic origin, correspondingly; Uij are the components of clastic deformation tensor, Cji are 
the coefficients of elasticity. It should be noted, that single-ion magneto-elastic interactions 
in Eq. (34) are written in the crystallographic coordinate systems XYZ, so, as distinct from 
Eq. (5), the indices of quadrupole moment components Q^'' = \{s^ s'' -\- s'^s^) are also defined 
in this system. 
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The values of elastic deformations, which appear because of spin configuration change, 
will be found by the minimaization of energies (33) and (34) by corresponding components 
of deformation tensor. As a result the following expressions are obtained: 



cii + ci2 - Zcfg/csa 



a/3 



a/3) ^ 

-°12 



''"(E^"/5«"«/3+E^3r'Qr (35) 



C33 



a/3 



En("/3) Z 
-°33 ■Sa'S/3J 

a/3 



a a/3 



X (^a - sis}) , (36) 



(fill + C12) 



C33 (Cll + C12) - 20^3 



Cl3 



. a/3 



33 *^<3f^ + E^33 ''a''/3 
a/3 



Cll + C12 



{2 ^ 7a/3SaS/3 + ^ (s^r^^ + B^^^) (Q^^ + Q^^) (37) 

a/3 a 

Q/3 

The Eq. (35) determines the isotropic striction of "easy" plane, or, what is the same, its 
expansion (or its contraction, depending on the signs of magneto-elastic coefficients), and 
Eq. (37) - the expansion/contraction along the crystal symmetry axis. 

The spontaneous deformation in singlet phase is defined by obtained Eqs. (35)-(37) after 
substitution in them corresponding values of spin variables: s = 0, = 0, Q^^ = Q^^ = 
1. After this it follows, that in this phase only the isotropic deformation of "easy" plane and 
expansion/constriction will be not equal to zero: 



''xx + '"yy ~ ~^ :~ 



B 



Cll +Ci2 - 2c?3/C33' 



,(0) 



Cl3 



B 



12 



(Cii + C12) 033 - 2cf3 ' 



(38) 



(39) 



where index (0) refers to the spontaneous magnetostriction. It can be seen, that in singlet 
phase the spontaneous deformations are defined only by single-ion magneto-elastic coefficients 
and satisfy the ratios: u^x = u^y = — 'uiz^C33/2ci3. The expressions (38) and (39) remain 
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correct in the magnetic field too, while h\\ < hg, i.e. in the region of the singlet phase 
stability. In other words, the striction, that is specified in this region by Eqs. (38)-(39), does 
not depend on the field. 

The induced striction appears only after the spin polarization occurrence, or in the fields 
/i|| > hs [47], and is described by Eqs. (35)-(37). They were written in the general form 
and included all admitted phenomenological parameters of magneto-elastic coupling, that 
has both exchange (inter-ionic) and single-ion (because of the change of crystal field, that 
affects the ions) origin. So, at the analysis of magnetostriction, it should be considered 
several, by our opinion, interesting cases. 

At first, let consider the magnetostiction, that is caused by isotopic exchange interaction. 
In the most of magnets, the corresponding magneto-clastic coupling docs not depend on spin 
directions and usually exceeds the anisotropic magneto-elastic one on more then order of 
magnitude. However, it is easy to convince, that despite the fact, that exchange magneto- 
elastic interaction does not depend on spin directions in crystal, the striction, that is 
generated by the external field, can be anisotropic. 

Indeed, it will be supposed, that in Eq. (34) only the magneto-elastic coefficients are 
finite Ai2 ^ 0, 711 ^ 0, and also C13 —^ 0. Such a situation can take place, for example, 
in the lamellar crystals. If magnetic sublattices are formed by spins in basal planes, then 
the intersublattice antiferromagnetic exchange depends basically on inter-atomic distances 
along the crystal symmetry axis. As for intrasublattice one, it depends on the inter-ionic 
distances in this plane. Then from Eqs. (38) and (39) one can find, that in singlet phase all 
Uj^j^ = and it is not infiuenced by the field. When spin polarization becomes finite, the 
exchange magnetostriction is represented by quite simple ratios: 

^ = s2 = s2 (40) 

*j^xx ^yy 

^=siS2 = s'cos2e (41) 

where 29 is, as above, the angle between sublattice spins, ul'-^^ = u^yP = — 2711/(011 -I- C12) 
and ull^f = — 2A12/C33 arc the values of induced striction at /ly = hfup. It should be noted, 
that for this case u^xx = u^yy — u^J = also. 

It follows from Eqs. (40) and (41), that in the region near hg of induced by external 

field phase transition there are singularities in striction field dependencies. The derivatives 
9uzz/dh\\ and duxx/dh\\ will have a jump at this point. 
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Phc. 5: Exchange striction, which is described by Eqs. (40) and (41), for parameters \ J\I D = 
0.05, I/D = 0.3, Jz/D = 1. The Une 1 corresponds to the "longitudinal" deformation, 
Uzz/ul''^^, and line 2 - to the "transversal" one u^xl'^xi^ at h || OZ. The line 3 corresponds 
to the "longitudinal" deformation at h _L OZ. 

The field behavior of induced exchange striction, that is defined by Eqs. (40) and (41), 
is shown on Fig. 5. It can be seen, that in the field h |j OZ the exchange striction (it is 
normalized and - depending on the sign of 711 - can be both positive or negative), that 
is caused by the intrasublattice interaction, does not change its sign and only "follows" 
the behavior s^(/i||). At the same time, the stiction, that is originated from intersublattice 
exchange (its sign depends on the sign of A12) for such a field orientation, is non-monotonic. 
This fact is a direct and simple consequence of change of cos 2^ sign. At first, the spin 
polarization grows (during the increasing of the field above /ig) practically at the opposite 
spin directions, so the striction (by absolute value) also increases. However, the further field 
growth gives rise to spin sublattice tilt, which becomes more and more noticeable, and in its 
own turn it causes the decreasing and passing through zero at ^ = 7r/4 of deformation. After 
such an angle configuration, when ^ — > 0, the magnetostriction again increases, reaching the 
maximum in the field /ly = /i/Hp- 

According to Eqs. (40) and (41) the field behavior of and u^z at h ± OZ will be 
similar. The exchange induced striction in this field changes monotonously, but begins to 
appear at the point /i^ = 0. It should be noted, that in the region /i^ ^ the striction 
is proportional to the Y?^. Here also the derivatives duzzjdh^ and duxxldh±_ changes 
continuously without jumps. 

Thus, the main distinction of induced striction at h || OZ and h _L OZ is, that in the 
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longitudinal field there should be a jump in the field derivative of striction behavior versus 
field, and in the transversal field this derivative changes continuously. 

Let note, that field behavior of magnetostriction, shown on Fig. 5, is qualitatively agreed 
with experimental data, which are obtained from the measurement of DTN deformation 
[22,23]. Indeed, in this compound the intersublattice exchange dominates and resulting from 
it the magneto-elastic couphng refers to chains Ni-Cl-Ni-Cl, which are parallel to the axis 
OZ. It is still unknown, whether Ni ions, which lie in the one basal plane, form the one 
sublattice, because the nearest chains are shifted on the half of a period along axis OZ. But, 
nevertheless, there are no doubts, that in this singlet magnet the determinant (together with 
single-ion anisotropy) is the intersublattice (antiferromagnetic) exchange and its anisotropy. 

It is also not excluded, that in singlet magnets the anisotropy of magneto-clastic interaction 
can be comparable with isotropic one. The anisotropic part of magneto-elastic coupling may 
include both exchange (interion) and single-ion parts [4]. Furthermore, it should be so in fact, 
if the magnetic system, to which DTN refers, is associated with strong single-ion anisotropy 
(which is the evidence of the essential spin-orbital interaction. In DTN, that is described 
by Hamiltonian (2), the single-ion anisotropy essentially exceeds the exchange interaction, 
the consequence of what is the formation of singlet ground state of the magnet in whole. 
So, in this case it is possible the next situation: the striction, that is caused by anisotropic 
interactions, exceeds the one, which is originated from isotropic exchange. 

As one more example, let consider such a case, when all magneto-elastic coefScients, 
except for B33 and B^^\ can be neglected. Then, if also C13 — > 0, the main will be, as it 
is seen from Eq. (27), the deformation of crystal along axis OZ: 
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This expression for striction can be written in the normalized form: 
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where according to the definition 
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is the striction in the field h\\ = hfup. 



(42) 



(43) 



It is shown on Fig. 6 the field behavior of induced striction, which is obtained, using 
Eq. (43). The most interesting case is h || OZ, when ^^3"*^ ^ 0, and ^^3^' -|- B^g^^ = 
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Phc. 6: "Longitudinal" Uzz/ulz^ magnetostriction versus field. Line 1 is satisfied at 
B^a"'^ ^ and B^"^ + B^f^ = 0, the line 2 - S^g"'' = 0> -^33^^ + ^33^^ 0, the line 3 - 
{Bg^' + <0 / (^33-^^ + <^ + Bg'^) = 2.5, Sr^V (4r^^ + <^ + = 
for h II OZ. The line 4 is obtained for the same parameter values as for line 3, but at h ± OZ 

0, to which refers the curve 1 on this figure. Corresponding magnetostriction is in direct 
proportion to the quadrupole moment Q^'^ . It should be said, that obtained magneto- 
elastic contribution is linear (but not quadratic, as it is usually) by magnetization, taking 
into account, that, quadrupole moment as function of field is similar (see Figs. 1 and 3) to 
the behavior of m(/i||). It should be also considered, that at the large exchange anisotropy 
Jz the magnetization versus — hs changes almost linearly, or m|| ~ — hg. So one can 
say, that the magnetostriction too (see line 1 on Fig. 6) behaves almost linearly in a whole 
region of antiforromagnctic phase existence. 

The field behavior of stiction at h 1| OZ, because of parameters choice, appears to be 
proportional to the m|, is corresponded with fine 2 on Fig. 6. 

There is shown also on Fig. 6 the example (fine 3) of common action of both anisotropic 
magneto-elastic striction mechanisms. At these ratios there takes place the competition 
between single-ion and inter-ion terms. In particular it is seen from fine 3 on Fig. 6, that the 
field behavior of striction, that is caused by the anisotropic magneto-elastic interactions, can 
be the same as the stiction, that is originated from the isotropic exchange, what is shown 
on Fig. 5. However, it follows form Eq. (43), that at h || OZ the single-ion and inter-ion 
contributions in striction can compete, and at h _L OZ there remains only the single-ion one. 
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The striction behavior in h _L OZ is shown as Hne 4 on Fig. 6, at the same parameters, at 
which the line 2 was obtained. Thus, in the considered examples of the competition between 
magneto-elastic interactions, it appeared, that in large fields h || OZ and h _L OZ the 
longitudinal striction has difTerent signs. 

In DTN the components of tensor of longitudinal striction have one sign and are close 
in the values [22]. Therefore, the assumption, that the induced longitudinal striction at the 
magnetization in fields h || OZ and h _L OZ , that are corresponded to h\\^±_ — > hfUp in this 
compound is caused by the intersublattice isotropic exchange interactions, is quite probable 
(believable). But, nevertheless, the additional investigations of magneto-elastic properties 
of the system are needed to prove unambiguously, that the observed deformation, under 
the effect of the field, is resulted from the intcrplanc (in DTN - intersublattice) exchange 
interaction only, and is not the the consequence of several field contributions, including from 
the spin quadrupole moment. 

7 Conclusions 

Thus, it was obtained, that the phase transition to magnetically ordered state, induced 
by magnetic field, in Van Vleck antiferromagnet, is the quantum phase transition. The 
spin polarization of the magnetic ion ground state is the order parameter of this phase 
transition, and for its description the Landau theory can be used. The considered transition 
is a consequence of competition of different interactions, and, what is important, it appears in 
the field, that is perpendicular to the easy plane. Such a field does not reduce the symmetry 
in this plane, leaving all directions in it equivalent. The conservation of degenaracy for 
sublattice magnetization directions in the easy plane is the crucial symmetrical condition of 
phase transition to the antiferromagnetic state with spontaneous magnetizations, that are 
lying in this plane. The account of stricition (see Eq. (36)) leads to the spontaneous lowering 
of the plane symmetry. 

It is also shown, that in the magnetic field induced antiferromagnetic phase the spin 
polarization (magnetization) of sublattice changes continuously from zero, reaching its maximum 
at the spin flipping point. As distinct from classical Neel antiferromagnets, in the magnetic 
ordered phase of Van Vleck (singlet) antiferromagnet the magnitude of sublattice magnetization 
strongly depends on the field. The same dependence on the field has an angle, that defines the 
deviation of sublattice magnetization from field direction. At the same time, the magnetization 
of a system as a whole, being weekly dependent from the field, shows almost hnear field 
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behavior (when the exchange anisotropy is also accounted). 

The calculations indicate, that in such an antiferromagnet the induced magnetostriction 
appears only in the magnetic phase. This magnetostriction in the small fields region is 
connected with the spontaneous formation of sublattice magnetizations. In the large fields 
(which are corresponded to spin fiiping field) the magnetostriction is basically determined 
by the sublattice magnetization rotation. 

Let emphasize two important issues. First one - is the possibility of induced striction, 
which is originated by the intrasublattice magneto-elastic interaction. In the classical antiferromagnets 
this part of induced magnetostriction of antiferromagnets is usually neglected, because of 
the paraprocess smallness. The second aspect - is connected with single-ion striction, the 
value of which is directly proportional to the spin quadrupolc moment; as a result, it occurs 
that the striction, which is caused by this quantity, has close to hnear dependence upon the 
field. 

Finally, the methodic notation should be made. The results, presented above, were 
obtained using the approximation of self-consistent field. It was supposed, that more accurate 
calculations will not bring any qualitative results, however they can infiuence quantitatively. 
Moreover, the magneto-elastic energy was written in the phenomenological form and contained 
a lot of parameters. So at the analysis of concrete compounds one should proceed from 
its characteristic hierarchy of interactions in the magneto-elastic energy, like it was made 
for example in this article. The separate paper will be devoted to detailed comparison of 
calculations with the available experimental data. 

We are grateful to Prof. S.M. Ryabchenko, who paid our attention on experimental 
articles [22,23] and the problem of magnetostricton properties of singlet magnets. 

CnHCOK jiHTepaTypbi 

[1] L. Neel, Magnetisme et champ molekulaire local, Conference Nobel, 11 decembere 1970. 
(Reimpression de Les prix Nobel en 1970, Stockholm) p. 57 

[2] A. S. Borovik-Romanov, Antiferromagnetism in Progress in Science [in Russian], (Izd- 
vo AN USSR, Moscow, 1962) 

[3] E.A.Turov, Physical Properties off Magnetically ordered Crystals (N.Y. -London, 
Academic Press, 1965) p. 222 



2fi 



[4] K.P.Belov, A.K.Zvezdin, A.M.Kadomtseva and R.Z.Levitin, Orientational Transitions 
in Rare-Earth Magnets (Science, Moscow, 1979) 

[5] A.F.Andreev, and V.I.Marchenko, Physics-Uspekhi, 23, 37(1980) 

[6] M.K.Wilkinson, J.W. Cable, E.O. Wollan, and W.C. Koehler, Phys.Rev., 113, 
497(1959). 

[7] M.E. Lines, Phys. Rev. 131, 546 (1963) 

[8] J.W. Leech, and A.J. Manuel, Proc. Phys. See, B59, 210 (1956). 

[9] C. Starr, F. Bitter, and A.R. Kaufmann, Phys. Rev., 58, 977 (1940). 

[10] D. Billerey, C. Terrier, A.J. Pointon, and J. P. Redoules, J.Magn.& Magn. Materials 21, 
187 (1980). 

[11] A.F. Lozenko, P.E. Parhomchuk, S.M. Ryabchenko, and P.A. Trotsenko, Low 
Temp.Phys., 14, 241 (1988) 

[12] V.M. Kalita, A.F. Lozenko, S.M. Ryabchenko, and P.A. Trotsenko, 
Jour.Exp.Theor.Phys, 99, 1054 (2004) 

[13] V.M. Kalita, A.F. Lozenko, S.M. Ryabchenko, and P.A. Trotsenko, Low Temp.Phys., 
31, 794 (2005) 

[14] B. Dorner, D. Visser, U. Stiegenberger, K. Kakurai, and M.Steiner, Z. Phys. B 72, 487 
(1988). 

[15] A. Harrison, and D. Visser, Condens. Matter 4, 6977 (1992) 

[16] S. Sachdev, Quantum phase transitions, Cambridge Univ. Press, Cambridge (1999) 

[17] R.S. Gekht, Physics-Uspekhi, 32, 233 (1989) 

[18] M.F. Collins, and O.A. Petrenko, Can. J. Phys., 75, 605 (1997). 

[19] B.S. Dumesh, Physics-Uspekhi, 43, 365 (2000) 

[20] V.S. Zapf, D. Zocco, B.R. Hansen, M. Jaime, N. Harrison, CD. Batista, M. Kenzelmann, 
C. Niedermayer, A. Lacerda, and A. Paduan-Filho, Phys Rev Lett., 96, 077204 (2006). 

27 



[21] S. A. Zvyagin, J. Wosnitza, C. D. Batista, M. Tsukamoto, N. Kawashima, J. Krzystek, 
V. S. Zapf, M. Jaime, N. F. Oliveira, Jr., and A. Paduan-Filho, Phys. Rev. Lett., 98, 
047205 (2007) 

[22] V.S. Zapf, V.F. Corea, CD. Batista, T. Murphy, E.D. Palm, M.Jaime, S. Tozer, A. 
Lacerda, and A. Paduan-Filho, J. Appl. Phys., 101, 09E106 (2007) 

[23] V.S. Zapf, V.F. Corea, P. Sengupta, CD. Batista, M. Tsukamoto, N. Kawashima, 
P. Egan, C. Pantea, A. Migliori, J.B. Bets, M. Jaime, and A. Paduan-Filho, E-print 
archives: preprint cond-mat/0705.0365vl (2007). 

[24] S.M.Stishov, Physics-Uspekhi, 47, 789 (2004) 

[25] Yu.V. Pereverzev and V.G. Borisenko, Sov. Phys. Solid State, 26, 1249 (1984) [in 
Russian] . 

[26] Ye.V. Rosenfeld, JETP-Letters, 24, 60 (1976) [in Russian] 

[27] V.M. Kalita, and V.M. Loktev, Jour.Exp.Theor.Phys, 98, 1006 (2004) 

[28] V.M. Kalita, and V.M. Loktev, Low Temp.Phys, 32, 114 (2006) 

[29] T. Haseda, N. Wada, M. Hata, and K. Amaya, Physica B, 108, 841 (1991). 

[30] Y. Tanaka, H. Tanaka, T. Ono, A. Oosawa, K. Morishita, K. lio, T. Kato, H.A. Katori, 
M.L Bartashevich, and T. Goto, E-print archives: preprint |cond-mat /0104287| (2001) 

[31] B.A. Ivanov, and R.S. Khymyn, Jour.Exp.Theor.Phys, 104, 307 (2007) 

[32] As pointed in the recent paper [31], the nematic (in the very meaning of this term) 
ordering in strongly anisotropic magnets is absent. However the ground state of such 
magnet, that is singlet, do not defer, in physical sense, from previous one, that is caused 
by the large biquadratic exchange [27,28]. 

[33] T. Nikuni, M. Osikawa, A. Oosawa, and H. Tanaka, Phys. Rev. Lett., 84, 5868 (2000). 

[34] M. Matsumoto, B. Normand, T.M. Rice, and M. Sigrist, Phys. Rev. Lett., 89, 077203 
(2003) 

[35] E. Ya. Sherman, P. Lemmens, B. Busse, A. Oosawa, and H. Tanaka, Phys. Rev. Lett. 
91, 057201 (2003). 



28 



[36] T. Radu, H. Wilheim, V. Yushanhai, D. Kovrizhin, R. Goldea, Z. Tylczynski, T. 
Luhmann, and F. Steglich, Phys. Rev. Lett., 95, 1272002 (2005). 

[37] A.I. Bugriy, and V.M. Loktev, Low Temp. Phys., 33, 37 (2007) 

[38] V.M. Loktev, and V.S. Ostrovskiy, Low Temp. Phys., 20, 775 (1994) 

[39] It should be denoted, that ground state functions, given in [20] for laboratory 
(crystallograpliic) coordinates system, can not be considered as eigen ones, until the 
self-consistent problem is solved (see Ref. [38]). 

[40] V.M. Kalita, I.M. Ivanova, and V.M. Loktev, Low Temp. Phys., 28, 475 (2002) 

[41] V.M. Kalita, and V.M. Loktev, Physics-Solid State, 45, 1523, (2003) 

[42] F.P. Onufrieva, JETP 89, 2270 (1985). [in Russian] 

[43] Yu. A. Fridman, Ph. Klevets, and O.V. Kozhemyako, J. Magn. Magn. Matter,. 264, 
111 (2003) 

[44] Yu.N. Mitsai, Yu.A. Fridman, O.V. Kozhemyako and O.A. Kosmachev, Low Temp. 
Phys. 25, 569 (1999) 

[45] E.A. Turov, and V.G. Shavrov, Physics-Uspekhi, 26, 411 (1983) 

[46] P.Morin, D. Schmitt, and E. du Tremolet de Lacheisseri, Phys. Rev. B21, 1742 (1980) 
[47] V.M. Kalita, I.M. Ivanova, and V.M. Loktev, Ukr. Phys. Jour. 50, 1159 (2005) 



29 



